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Abstract We study the polarization efficiency (defined as the ratio of polarization to extinc¬ 
tion) of stars in the background of the small, nearly spherical and isolated Bok globule CB4 
to understand the grain alignment process. A decrease in polarization efficiency with an in¬ 
crease in visual extinction is noticed. This suggests that the observed polarization in lines of 
sight which intercept a Bok globule tends to show dominance of dust grains in the outer lay¬ 
ers of the globule. This finding is consistent with the results obtained for other clouds in the 
past. We determined the distance to the cloud CB4 using near-infrared photometry (2MASS 
JHKS colors) of moderately obscured stars located at the periphery of the cloud. From the 
extinction-distance plot, the distance to this cloud is estimated to be (459 ± 85) pc. 
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1 INTRODUCTION 

Linear polarization measurement of background field stars in a star forming cloud has been used to study 
the large scale structure of the interstellar magnetic field (Heiles (2000); Fosalba et al. (2002); Bertrang 
et al. (2014)). The efficiency with which polarization is produced by dust grains is sensitive to physical 
conditions of the cloud, with a tendency to drop systematically with increasing extinction or optical depth 
inside interstellar clouds (Goodman et al. (1995); Gerakines et al. (1995); Whittet et al. (2008); Chakraborty 
et al. (2014)). This observation suggests that the observed polarization in a particular direction or line of 
sight which intercept a dark cloud show dominance of dust grains in the outer layers of the cloud. However 
background stars observed through dense region of the cloud show smaller degree of polarization which 
suggests much lower dust columns in the observed lines of sight. This also suggests that the polarization 
studies of background stars in a cloud can well explain the external magnetic field, but it can not provide 
information on magnetic field structure within the cloud (Whittet (2005)). 

To model the observed variation in polarization with extinction along different lines of sight, Jones 
(1989) and Jones et al. (1992) studied this in terms of a magnetic field with distinct random and uniform 
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increase in extinction (Tamura et al. (1987); Goodman et al. (1992); Gerakines et al. (1995); Whittet et 
al. (2008); Chakraborty et al. (2014)). This can be represented by a power law p/Ay oc Ay° for some 
power law index a. From the study of dust grain alignment in Taurus dark cloud, Gerakines et al. (1995) 
found a strong correlation between polarization efficiency and extinction, which was fitted by power law 
p/Ay oc Ay 0 56 , in the magnitude region (0 < Ay < 25). Whittet et al. (2008) also studied the polarization 
efficiency toward molecular clouds with an aim to understand grain alignment mechanisms in dense regions 
of the interstellar medium and found an empirical relation: p\/r\ oc [Ay ]~ 0 - 52 (where t\ is the optical 
depth). They did not find any significant change in the behavior in the transition region between the diffuse 
outer layers and dense inner regions of clouds. Recently, Chakraborty et al. (2014) studied the polarization 
efficiencies in three selected Bok globules CB56, CB60 and CB69. They found a trend of decreasing polar¬ 
ization efficiency with increasing extinction which can be well represented by a power law, py/Ay oc Ay a , 
where cr = -0.56 ± 0.36, -0.59 ± 0.51 and -0.52 ± 0.49 for CB56, CB60 and CB69, respectively. To 
study the alignment of grains, radiative torque mechanism has been used by researchers to calculate the ex¬ 
pected relationship between polarization efficiency and extinction through a homogeneous cloud (Dolginov 
& Mytrophanov (1976); Draine & Weingartner (1996); Draine & Weingartner (1997); Lazarian & Hoang, 
(2007); Whittet et al. (2008); Hoang & Lazarian (2014)). Hoang & Lazarian (2014) studied grain alignment 
by radiative torques for various environment conditions, including interstellar medium, dense molecular 
clouds and accretion discs. 

The determination of distance to an interstellar cloud help researchers to obtain luminosities of the 
protostars embedded in the clouds (Yun & Clemens (1990)) and to estimate the sizes, masses and densities 
of the cloud (Clemens et al. (1991)). Several methods have been adopted to estimate the distance to cloud. 
These are Wolf diagrams method (Wolf (1923)), star counts method (Bok & Bok (1941)), spectroscopic 
method (Hobbs et al. (1986); Hearty et al. (2000)), photometric method (Snell (1981); Peterson & Clemens 
(1998); Maheswar & Bhatt (2006); Maheswar et al. (2004); Maheswar et al. (2010)) etc. Maheswar & Bhatt 
(2006) determined the distances to nine dark globules by a method which uses optical and near-infrared 
photometry of stars projected towards the field containing the globules. Recently, Eswaraiah et al. (2014) 
determined the distance to starless dark cloud LDN 1570 using near-IR photometry from 2MASS. 

In this work, we plan to study the polarization efficiency of field stars background to the isolated Bok 
globule CB4. We also determined the distance to the cloud CB4 using a technique developed by Maheswar 
et al. (2010) which uses JHK near-infrared photometry. 

2 OBJECT 
2.1 CB4 

CB4 is a small, nearly spherical and isolated Bok globule situated at a distance of 600 parsec (Dickman & 
Clemens (1983)). CB4 is situated far away from the galactic plane (b « —10°). The angular size of this 
cloud on the POSS maps is 2' x 1' and the position angle of the major axis (measured east from north) is 45°. 
An IRAS point source is located at southern edge of the cloud which is detected at 100 pm only (Clemens 
& Barvainis (1988)). No further evidence of an embedded source has been observed up to 160 pm. CB4 
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(1991)). The observations have shown that the globule velocity is about —11.3 km s -1 with the velocity 
of dispersion 0.57 km s -1 . Furthermore, the study of deep IRAS aperture photometry and molecular line 
showed that CB4 is cold (« 7K) and does not contain a dense core (n > 10 4 cm -3 ) (Clemens et al. (1991), 
Kane et al. (1994) and Turner et al. (1997)). 

3 DATA 

3.1 Polarimetric data 

The polarimetric data of 80 field stars of CB4 at V-band has been taken from Kane et al. (1995). Four fields 
were initially observed by them to cover a region of ~ 8' x 8' surrounding CB4. Two other fields were also 
observed again with a slightly larger pixel size and using longer integrations. The histogram of polarization 
position angles (PA) of 80 field stars displayed a strong peak around mean PA 65° with a dispersion of 
13°. This position angle represents the preferred magnetic field direction and the narrow width of the peak 
indicates a larger uniform field. The unweighted mean polarization and weighted mean polarization were 
calculated by Kane et al. to be (2.84±0.25)% and (0.96±0.07)% which indicate a large uniform magnetic 
field. 

3.2 Photometric data 

The near IR photometric data (J, H and K $) of the field stars of CB4 has been taken from the 2MASS 
All-Sky Catalog of Point Sources (Cutri et al. (2003)) which satisfies the following criteria: 

1. photometric uncertainty aj, an, &k s < 0.035 in all three filters, 

2. signal-to-noise ratio (SNR) > 10 which is shown by photometric quality flag of “AAA” in all three 
filters. 

The position of 80 stars are shown by RA(2000) and DEC(2000) shown in Table 2 of Kane et al. (1995). 
It is to be noted that we will designate all field stars by its serial number used in Kane et al. (1995)’s paper 
and this numbering has been shown in column number two of Table 1. The photometric magnitude of each 
field stars has been searched via CDS’s VizieR service, using 2 arcmin search radii as query. The J, H 
and Ks magnitudes of 80 stars are available in the catalogue. We also checked the SIMBAD astronomical 
database to search for any variable star out of 80 field stars, but we could not find any variable star. We found 
that only 60 stars satisfy the above criteria and our analysis is restricted to these stars only. Further we could 
not find the uncertainties of magnitudes of six stars in the 2MASS catalogue which are star number 4, 20, 
27, 50, 51 and 74. So we discarded them from our analysis. Thus total number of stars considered for our 
analysis is 54. 

4 RESULTS 

4.1 Spectral types and Ay 


We used the technique developed by Maheswar et al. (2010) to infer spectral type, and hence absolute 


4 


Barman & Das 


gives distances to globules that are relatively close (< 500 pc) with a precision of ~ 20%. In this technique, 
the fields containing the cores are divided into small sub-fields to avoid complications created by the wrong 
classifications of giants into dwarfs. Actually the increase in the extinction caused by a cloud should occur 
almost at the same distance in all the fields, if the whole cloud were located at the same distance, the 
wrongly classified stars in the sub-fields would show high extinction not at the same but at random distances 
(Maheswar et al. (2010)). But for small and isolated clouds, it would be difficult to divide the field containing 
the cloud into sub-fields with a sufficient number of stars projected on to them. Therefore the distance 
measurement is somehow hard for CB4 as it is a small, nearly spherical and isolated Bok globule. It would 
be difficult to subdivide the field with only limited number of stars. So we consider a single field surrounding 
the core of the CB4. 

In this technique, stars with (J — Kg) < 0.75 were selected from the fields containing the cloud to 
eliminate M-type stars from the analysis. Then a set of dereddened colors ( J — H) 0 and ( H — Kg) a has 
been produced for each star from their observed colors (J — H) and ( H — Kg) using trial values of Ay 
and a normal interstellar extinction law (i.e., total-to-selective extinction value, Ry = 3.1) in the equations 
(Rieke & Lebofsky (1985)): 

(J — H) a = (J - H) - 0.107 x A v (1) 

{H - K s )o = {H - K s ) - 0.063 x Ay (2) 

The maximum extinction that could be computed using this method is limited to Ay ps 4 mag. The 
trial values of Ay is chosen in the range 0 — 10 mag with a step size of 0.01 mag. The calculated set of 
dereddened color indices are then compared with the intrinsic color indices ( J — H)i n and (II — Kg)i n of 
normal main sequence stars and giants, produced using the procedures discussed in Maheswar et al. (2010). 
The best fit values of the dereddened colors to the intrinsic colors giving a minimum value of x 2 then give 
the corresponding spectral type and Ay for the star. The solutions which give Xmin — 0-1 is considered for 
this analysis. 

The uncertainty in Ay is given by 

cr(Ay) = ^A.7 2 .a 2 JH + 7.9 2 .a 2 HKs + 2 x 37 x cov{JH , HKg) (3) 

where a 2 H = a 2 + <j 2 h , crj IKs = + a\ s and cov(JH, HKg) = r s x aju^HKs ■ The Spearman rank- 

order correlation coefficient (r s ) is calculated from uncertainties in (J — H) and (H — Kg) colors which 
shows a strong correlation between them. The maximum uncertainty in the Ay in our case is estimated to 
be 0.5 mag. 

We found that 14 stars have ( J — Kg) > 0.75, so we also discarded them from our analysis. We used 
this technique to estimate the Ay and absolute magnitudes for 40 stars. It has been found that only 20 stars 
have Xmin — 0-1 f° r which it is possible to infer the spectral type and the Ay. The results are shown in 
Table-1. 

In Maheswar et al. (2010)’s technique, the field stars in the globule are assumed to be normal main 
sequence stars and therefore the luminosity class V has been assigned to all stars. However the presence 
of non main sequence stars may result wrong distance to stars. The intrinsic color and luminosity may be 
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plot ( B — V) vs ( V — R ) color - color diagram. The photometric data ( BVR ) of the field stars of CB4 has 
been taken from the VizieR database of astronomical catalogues, namely, UCAC4 (Zacharias et al. (2013)) 
and NOMAD (Zacharias et al. (2005)). The B 1 V. R data of star with serial number 79 is not available in 
the catalogue. In Fig. 1, color - color diagram along with error bars for 19 field stars of CB4 has been 
plotted. It can be seen from the graph that most of the field stars are clustered together along a diagonal 
locus. These diagonally clustered stars in the colour - colour diagram represent the main-sequence stars, 
while the scattered ones (Star # 32, 66 and 76) represent non-main-sequence stars. So we discarded three 
stars from our analysis. To calculate r s , we plotted uncertainties in (J — H) and (II — Ks) colors for 16 
stars which is shown in Fig.2. The value is found to be 0.95. The results obtained from this technique are 
shown in Table-2. 
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Fig. 1 Plot of B — V versus V—R colour - colour diagram of 1 9 field stars in the CB4. Diagonally 
clustered points represent stars that are in the main sequence. The stars that do not belong to the 
main sequence are represented by their serial numbers as mentioned in Table 1. 


4.2 Polarization efficiency 

The study of the variation of background star polarization with extinction (Ay) and polarization efficiency 
(defined as the ratio of p to Ay) can provide useful information regarding the nature of dust and the magnetic 
field associated with the cloud (Gerakines et al. (1995); Arce et al. (1998)). The polarization efficiency 
depends on both the properties of the dust grains and the degree of alignment of the grains. Whittet et 
al. (2003) calculated the theoretical upper limit which is given by p/A v < 14% mag -1 . However, the 
observational upper limit of p/A v is found to be ~ 3% mag -1 , a factor of four less than the predicted value 
for the ideal scenario. 

In this work, we will study the polarization efficiency of background field stars of CB4 to investigate the 
alignment of non spherical dust grains. It is to be noted that we designated all field stars by its serial number 
used in Kane et al. (1995)’s paper and this numbering has been shown in column number two of Table 1. 










6 


Barman & Das 


0.04 


0.035 - 


X 

o’ 


0.03 


0.025 


r., = 0.95 


0.025 


0.03 0.035 0.04 0.045 


HKo 


Fig. 2 The relation between the cthk s and 17 J n of the 16 stars. The r s value indicates a strong 
correlation between the uncertainties of the two colors. 


Table 1 The J, H and K$ magnitudes of 20 stars along with spectral type, absolute magnitude 
(Mk s ) and extinction (A,,). The second column refers to the star number which corresponds to 
the serial number in table 2 of Kane et al. (1995). 


S/N 

Star 

J 

(mag) 

H 

(mag) 

I<s 

(mag) 

Sp Type 

M Ks 

A v 

1 

1 

10.942±0.021 

10.604i0.018 

10.533i0.018 

G4V 

3.45 

0 

2 

9 

11.026±0.021 

10.461i0.019 

10.305i0.018 

K4V 

4.22 

0.40 

3 

11 

11.888±0.021 

11.412i0.019 

11.282i0.018 

K2V 

4.15 

0.23 

4 

14 

11.555±0.021 

11.185i0.018 

11.099i0.018 

G9V 

4.03 

0 

5 

17 

11.653±0.021 

11.186i0.019 

ll.li0.018 

K1V 

4.13 

0 

6 

18 

12.715±0.022 

12.515i0.023 

12.396i0.023 

A8V 

1.64 

1.05 

7 

23 

13.083±0.021 

12.702i0.023 

12.607i0.025 

K0V 

4.15 

1.34 

8 

26 

11.52±0.021 

11.206i0.019 

11.108i0.018 

G3V 

3.34 

0.24 

9 

32 

12.978±0.022 

12.552i0.023 

12.382i0.023 

G2V 

3.29 

1.44 

10 

35 

10.335±0.021 

10.199i0.021 

10.127i0.018 

F0V 

1.86 

0.16 

11 

36 

13.332i0.022 

13.036i0.029 

12.98i0.031 

G2V 

3.29 

0 

12 

41 

12.897i0.026 

12.354i0.029 

12.268i0.029 

K3V 

4.15 

0 

13 

52 

11.055±0.021 

10.863i0.019 

10.773i0.02 

F1V 

2.00 

0.45 

14 

61 

12.691i0.021 

12.418i0.021 

12.321i0.025 

F7V 

2.75 

0.44 

15 

62 

13.038i0.025 

12.693i0.027 

12.566i0.029 

GOV 

3.19 

0.81 

16 

66 

13.19i0.025 

12.865i0.031 

12.789i0.034 

G4V 

3.45 

0 

17 

71 

12.091i0.021 

11.881i0.021 

11.783i0.023 

F1V 

2.00 

0.59 

18 

76 

12.551i0.022 

12.238i0.023 

12.104i0.024 

F6V 

2.62 

0.98 

19 

77 

13.3i0.024 

12.9i0.027 

12.784i0.03 

G9V 

4.03 

0.35 

20 

79 

13.483i0.026 

13.103i0.026 

13.037i0.032 

G9V 

4.03 

0 


plotted Ay versus p for all the field stars in the cloud CB4. The solid line in the figure represents optimum 
alignment in the general interstellar medium according to Serkowski et al. (1975): p max /A y = 3 percent 
mag -1 , where Ry = 3.1. It can be seen that only four stars located outside the limit and majority of stars 
lie below the solid line which suggests that the characteristics of material composing CB4 is consistent with 
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Table 2 Polarization ( p ), extinction (Ay), polarization efficiency ( p/Ay ) and distance (in parsec) 
of 16 background stars in CB4. The second column refers to the star number which corresponds 
to the serial number in table 2 of Kane et al. (1995). 


S/N 

Star 

P 

(%) 

Ay 

p/A v 

(% mag -1 ) 

d 

(pc) 

1 

1 

0.30±0.05 

0 

- 

261 ±48 

2 

9 

1.66±0.41 

0.40 

4.15±0.53 

161±30 

3 

11 

1.26±0.38 

0.23 

5.48±0.51 

264±49 

4 

14 

0.49±0.08 

0 

- 

259±48 

5 

17 

2.85±0.38 

0 

- 

248±46 

6 

18 

1.52±0.06 

1.05 

1.45±0.41 

1342±249 

7 

23 

1.37±0.12 

1.34 

1.02±0.43 

459±85 

8 

26 

1.07 ±0.04 

0.24 

4.46±0.34 

353±65 

9 

35 

0.27±0.07 

0.16 

1.69±0.36 

446±83 

10 

36 

1.16±0.15 

0 

- 

867±162 

11 

41 

1.94±0.23 

0 

- 

420±78 

12 

52 

0.97±0.12 

0.45 

2.16±0.37 

555±103 

13 

61 

0.47±0.08 

0.44 

1.07±0.40 

802±149 

14 

62 

1.14±0.27 

0.81 

1.41 ±0.55 

720±134 

15 

71 

0.55±0.12 

0.59 

0.93±0.40 

878±163 

16 

77 

1.51±0.25 

0.35 

4.31±0.54 

553±103 


In Fig. 3b, we plotted Ay vs p/Ay for cloud CB4. The plot shows that the polarization drops with 
increasing extinction. Thus background stars observed through dense region of the cloud show smaller 
degree of polarization suggesting much lower dust columns in the observed lines of sight. This suggests 
that the observed polarization at optical region in a particular line of sight which intercept a dark cloud 
show dominance of dust grains in the outer layers of the cloud. Recently, Chakraborty et al. (2014) found 
a trend of decreasing polarization efficiency with increasing extinction which can be well represented by a 
power law, py /Ay oc Ay a , where a = —0.56 ± 0.36, —0.59 ± 0.51 and —0.52 ± 0.49 for Bok globules 
CB56, CB60 and CB69, respectively. It is to be noted that the study of polarization efficiency at optical 
region can give us information about the grain alignment in outer layers of the cloud. The alignment of 
grains by radiative torque has been used by researchers to calculate the expected relationship between 
polarization efficiency and extinction through a homogeneous cloud (Dolginov & Mytrophanov (1976); 
Draine & Weingartner (1996); Draine & Weingartner (1997); Lazarian & Hoang, (2007); Whittet et al. 
(2008); Hoang & Lazarian (2014)). 

4.3 Estimation of distance 

The photometric distance d to a star is estimated using equation: 

d (in pc) = iq( k s-m Ks+ 5-a K s )/5 (4) 

where Ks, Mk s and Ak s are apparent magnitude, absolute magnitude and extinction, respectively. The 
distance to the cloud is typically estimated from the first star that shows a sharp peak in extinction in an Ay 
vs. d plot. It is to be noted that Ak s = 0.112 x Ay. 

The uncertainty in distance is estimated using the expression (Maheswar et al. (2010)), 


n 
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Fig. 3 (a) Plot of Ay versus p for 11 field stars. The solid line in the figure represents optimum 
alignment in the general interstellar medium according to Serkowski et al. (1975): Pmax/Ay = 3 
percent mag -1 , where Ry = 3.1. It can be seen that four stars located outside the limit, (b) plot 
of Ay versus p/Ay (where Ry = 3.1) for 11 background stars (see Table 2). 



Fig. 4 Extinction (Ay) versus distance for stars in the vicinity of the CB4. The average extinction 
(< Ay >) of all the field stars is 0.36 and is shown by solid line. The vertical dashed line is drawn 
at a distance of 459 parsec where sudden rise in the Ay occurs. 

where <tk s is the uncertainty in Ks band, <Jm Ks is the uncertainty in the estimation of the absolute mag¬ 
nitude and <Ja Ks is the uncertainty in the Ak s estimated by the method. The uncertainty in the <Jm Ks is 
assumed to be 0.4 while calculating er^ in the distances for all the stars. The typical errors in our distance 
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The extinction (Ay) versus distance ( d ) plot has been extensively used by several investigators 
(Dickman & Clemens (1983); Kenyon et al. (1994); Peterson & Clemens (1998); Maheswar et al. (2004); 
Maheswar et al. (2010) etc.) to determine the distance of globules. It is noticed that the extinction of stars 
located in front of the dark globule shows almost zero reddening value. The presence of cloud shows sudden 
rise of extinction at a certain distance. The distance to the cloud is typically estimated from the first star 
that shows a sharp peak in Ay. We notice from Fig. 4 that star number 23 shows a sudden rise of extinction 
and then a fall in extinction is observed. Thus the distance of the cloud is (459 ± 85) parsec. Using the 
foreground star method and color excess approaches, Dickman & Clemens (1983) estimated distance to 
the CB4 to be ~ 500 — 700pc. They therefore adopted 600 pc as the distance to this globule. The distance 
determination using JHK near-infrared photometry is slightly smaller than the distance derived by Dickman 
& Clemens (1983). However, the present analysis is restricted to a limited number of data points. 

5 SUMMARY 

1. The polarization efficiency of CB4 is studied and shows a decrease in polarization efficiency with an 
increase in extinction along the observed line of sight. This finding is consistent with the results obtained 
for other clouds in the past. 

2. The method for determining distances to Bok globules using the NIR photometric technique has been 
described. The plot of Ay versus d has been made and the distance at which Ay shows a sudden rise 
of extinction is taken as the distance to the cloud. From the study of the distribution of reddening, we 
estimate a value of (459 ± 85) pc for the distance to CB4. 
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